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Using exact Lanczos diagonalizations we have shown that 
the pair-hopping model for negative W exhibits a phase tran- 
sition into n- superconducting state. The transition occurs 
at any band filling and the critical value W c varies between 
W c — — 1.75i at half-filling to W c — —It for two particles on 
the lattice. 

PACS numbers: 71.27. +a. 71.28.+d, 74.20.-z 



There is considerable interest in the pair-hopping 
model as a phenomenological model which captures some 
of the essential physics of the high-T c superconducting 
materials. The model was proposed ten years ago by 
Penson and KolbEl as a model exhibiting the real-space 
pair formation mechanism. The Hamiltonian of the pair- 
hopping model contains, in addition to the usual one- 
electron hopping term, a term that hops singlet pairs of 
electrons from site to site and in the one-dimensional 
case is given by 

ti = ~ t ^ ~ ' J { C n,a Cn +' s -, a + C n+l,a C n,a) 
n,ct 
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(c n ,a) is the creation (annihilation) operator 
for an electron with spin a at site n. Thus t and W are 
the single-electron and pair-hopping amplitudes respec- 
tively, and the competition between these two sources 
of dekjcalization leads to a rich phase diagram of the 
modeltTu. For W = 0, we have a system of free elec- 
trons, whose properties are exactly known. In the oppo- 
site limit t = all sites are doubly occupied or empty, the 
model is superconducting by construction and equivalent 
to the XY modelEl. The correspondance with the XY 
model can be seen by introducing S — 1/2 pseudospins 
T n , with T+ = 4 )T 4 4 ,T* = (4 T c„ jT +cj u c, a - l)/2. 
Then the Hamiltonian becomes 



n = -wY,(T+T n 



h.c). 



(2) 



This picture holds true for either sign of W, but it is im- 
portant to note that W — > — W is not a symmetry of the 
model ([!]) and the way the system reaches its limiting 
behaviour at \W\ » t is genuinely different for nega- 
tive and positive W casescl. Moreover, as we show in this 
paper, if for positive— W a singlet superconducting state 



with order parameter Ass = ( c 



t t 



corresponding 



to the usual Cooper pairs is realizedcL in the case of 
negative-W , an i]- superconducting stateO with order pa- 
rameter = (cj. i-cj._ fc I ) corresponding to the pairing 
with total momentum equal to n is the ground state of 
the system at W < W c w -1.75*. 

The model was mainly studied in the case of half-filled 
band and W > 0. Penson and Kolb used exact diago- 
nalizations data for chains up to 10 sites and found a 
phase transition into a superconducting state at which 
the spin gap (or single-particle excitation gapl opens for 
W > W c ~ 1.4*. Later Affleck and MarstonH analysed 
the model within the framework of the weakly-coupling 
continuum limit approach. They found that for any 
W > the spin excitation spectrum is gapped, while 
the charge excitation spectrum is gapless and the singlet 
superconducting instabilities are most divergent in the 
ground-state. In the case of W < (| W |<SC *) at half- 
filling there is a gap in the charge excitation spectrum, 
the spin sector is gapless and the ground state corre- 
sponds to an insulator. Moreover Affleck and Marston 
argued the absence of any other transitions for W > 
and the necessarity of a phase transition into a sector 
with gapped spin excitations and gapless charge degrees 
of freedom for — W ^> t. This scenario of the ground- 
state phase-diagram was recently confirmed by Sikkema 
and Afflecko who used the Density Matrix Renormaliza- 
tion Group technique (DMRG) for chains up to 60 sites. 
For W < 0, they found a transition into a spin gapped 
phase at W < W c w -1.5*. 

In this paper we focus our attention on the transition 
into an ^-superconducting state for negative-W . This 
transition corresponds to a drastic change of the ground- 
state, after which the one particle hopping term is almost 
frozen out. To prove this, we used exact Lanczos diag- 
onalizations. We have shown that the transition into r\- 
superconducting state takes place at arbitrary filling and 
weakly depends on the band-filling. Our data suggests 
that the finite size effect are extremely small, that is why 
we restricted ourself to systems up to 10 sites. 

As far as the main phenomenon characterizing the 
transition into an ^-superconducting state takes place 
in the momentun space it is convenient to rewrite the 
Hamiltonian in this space 

H = -2tJ2 A, a c k, a cos(fc) - 2W J2 a q a Q C ° S (Q) (3) 

k,a Q 

where Aq — c\ |Cg_ fc ^ is the creation operator of 

pair of electrons with opposite spins and total momentum 
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Q, L is the size of the system. As far as the total momen- 
tum of the system Qtot is conserved we can consider each 
sector of the Hilbert space for a given Q to t independently. 
It is clear that the ground state of the system belongs to 
the sectors Qtot = or Q to t — t. States with Qtot 7^ 
and tt exhibit a broken time reversibility symmetry and 
are excited states. _ 

As it was shown by Yangp the general condition ensur- 
ing the possibility for 77-pairing is 

e(fc) + (.{tt — k) — constant (4) 

where e[k) is the bare electron spectrum and tt — 
(tt, ..,7r). The particularity of the pair-hopping model 
could be observed even in the case of two particles on 
a lattice. Below we consider the one-dimensional case, 
however the arguments are valid in any dimensions if the 
bare electron spectrum ensures the condition (^) . 

Let us first analyse the W < case. For W = the 
ground-state energy is Eq = —At and the total momen- 
tum is Qtot — 0. As we switch on W, it can be shown eas- 
ily that the lowest energy in the Qtot = subspace goes 
to -2W > (for \W\ > t). On the other hand, the low- 
est energy in the Qtot = tt sector is 2W < 0. Thus there 
should be a transition at some critical value, from the un- 
paired state in the Qtot = sector into the 77-pairing state 
in the Qtot — tt sector. We can roughly estimate the criti- 
cal value to be W c ~ — It. Indeed after the transition the 
ground-state corresponds to the wave function A\\Q >, 
which clearly describes an 77-pair. It is clear that after 
the transition the one particle hopping term is completely 
frozen out (i.e < n{k) >= J2 a < c \ a c ka >— constante) . 

In the opposite case W > 0, the ground-state always 
remains in the Qtot — subspace and its energy continu- 
ously goes from -At at W = to -2 FT for W > t. There 
is no more transition for W > and the weight of Cooper 
pair continuously increases up to 1 when W = +00. In 
this limit the ground state wave function is ^4q |0 >. Note 
that only for W = ±00 the ground-states are equivalent 
up to a trivial 7r/2 shift of all momenta, reflecting the 
W — ► —W symmetry of the XY model. 

As we show below using exact diagonalizations data 
the picture essentially remains in the case of many par- 
ticle systems. Namely for W < 0, there is a phase 
transition (presumably of the first order) at the critical 
value of the pair-hopping amplitude from the insulating 
(1/2-filled case) or Luttinger-liquid state (lower filling) 
into a ^-superconducting state. The transition occurs 
at any band filling and the critical value W c varies be- 
tween W c — — 1.754 at half-filling to W c ~ — It for two 
particles on the lattice. The transition to 77-pairing is 
characterized by the strong reduction of the one particle 
hopping term. However, it should be stressed that due 
to the quasi-bosonic character of the pairs, the weight of 
unpaired particles is extremely small but finite after the 
transition. 

In the following part, we will present numerical data 
considering the lowest energy in the Q = and Q = tt 



sectors. In Fig.l we have plotted the lowest energy in 
each sectors at half- filling, as a function of W for two dif- 
ferent cases L — 8 and L — 10. In Fig. la we clearly ob- 
serve a different behaviour of the ground-state energy for 
W > and W < 0. In the case W > the ground state 
energy continuously goes to the XY model ground-state 
energy (| W \— 00), and the ground-state always remains 
in the Q = sector. There is no trace of any additionnal 
transition for W >.D, this is in complete agreement with 
the DMRG resultsu. Let us now consider the W < 
case. Below a critical value W c ~ —1.754 the energy is 
quasi linear in W and becomes very close to the ground 
state energy of the XY model. Hence after the transition 
the one particle hopping term is almost frozen out. This 
transition can be observed more clearly in Fig. lb. Indeed 
before the transition the total momentum of the ground 
state is and is tt after the transition. This transition 
from the Q = to Q = tt subspace depends only on the 
parity of the number of pairs, therefore the transition is 
easily observed when this number is odd. 

To distinguish superconducting phases corresponding 
respectively to IT > and W < 0, we calculated the 
distributions of < > and < AqAq > in both cases 
(Fig. 2). For W > the < > distribution continuously 
goes to the limiting case < >— 1 (i.e. \W\ = 00) 
and the distribution of < AqAq > has a strong peak at 
Q = 0. In the opposite case, for W < W c , < nk >— 
1 (reflecting the frozening out of one electron hopping 
term) and < AqAq > shows a peak at Q — tt. For 
W > 0, pairs with Q = appear continuously in the 
system. On the other hand, when W < the Q = tt pairs 
appear spontaneously when W < W c . The competion 
between the i-term and the FT-term leads to the total 
destruction of the old band structure at the critical point. 

To show that Aq with Q = or tt is a proper operator 
to describe the physics of the system we calculate, 

z _ |<iV e -2|A Q |Y e >| 
Q {<N e \A Q A Q \N e >y/* 

where \N e > is the ground-state with N e particles. This 
overlapp measures the weight of quasi-particle with total 
momentum Q in the ground-stateo. We checked numer- 
ically that Zq is zero for all values of Q except and 
tt. In Fig. 3 we calculated this quantity for Q = and 
Q — tt. For W > 0, Z n — whilst Zq is finite and 
continuously approaches its limiting value. The jump at 
W = is a finite size effect, it can be shown that it goes 
as 1 / \[Tj. Thus the usual Cooper pairs appear continu- 
ously for IT > 0. 

In the opposite case W < 0, Z = 0. Depending on the 
parity of the number of sites Z„ is or 1 for W = 0~. 

For L = 6, Z^ = before W = —1.54, the jump which 
appears at this point corresponds to a 'crossing' with 
some excited state, while the final transition correspond- 
ing to ^-pairing occurs only at W ~ —1.754, after this 
transition Z^ ~ 1. 
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In the case of L = 8, Z v = 1 for W — 0~ this is 
common when L = AL Q due to the presence of electrons 
with momentum pp — ir/2. When we increase \W\ the 
probability to find a particle on the Fermi surface is re- 
duced, this is reflected by the reduction of Z„. After the 
transition to 77-pairing ~ 1 . This suggests that the n— 
pairs appear spontaneously in the ground-state after the 
transition. And Aq and A v are proper operators to de- 
scribe the superconducting states for W > and W < 
respectively. 

Another way to visualize this transition is to show the 
pairing phenomenon in real space. The correponding 
quantity is given by, 



1 L 

z=l 



< nil >< n il > 



(6) 



it is plotted in Fig. 4a. 

This picture shows that for W > the pairs appear 
continuously in the system. In the opposite case W < 0, 
and \W\ « f the tendancy to pairing is reduced, this 
is in agreement with the effective repulsive character of 
the pair-hopping interaction in this limitcl. This effective 
repulsive character remains up to \ W\ < 0.54 after which 
a tendancy to pairing appears. However the transition to 
ry-state corresponds to the jump in F(W) at W ~ — 1.75i. 
As it is clear from this figure that the finite size effects 
have only a weak influence on this value. 

To show that the physics does not depend on the band 
filling, we analyse the model away from half-filling. In 
Fig. 5 we plotted the lowest energy as a function of W in 
two particular cases. It is clear from this picture that the 
transition to restate remains away from half-filling. To 
analyse the band-filling dependance of the critical value 
we calculated F(W) for a given size L = 10 and differ- 
ent number of particles (fig4b). This picture suggests 
that the critical value weakly depend on the band filling. 
When we increase the band filling W c goes from -2t (two 
particles) to —1.754 (half- filled case). 

To conclude, in this paper we have shown that the 
pair-hopping model for W < exhibits a phase transi- 
tion into the ^-superconducting state. The critical value 
at which the transition takes place weakly depends on 
the size of the system and on the band-filling, it varies 
between —24 < W c < —1.754. For W. > our results are 
in agreement with previous worksau, and shows a con- 
tinuous second-order transition to usual superconducting 
state at W = + , with no additionnal transition for any 
W > 0. We argue that this phenomenon will remain un- 
changed in higher dimensions. Investigation of the 2D 
case is under consideration. 
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FIG. 1. Lowest energy vs. W calculated in the Q — 
(open circles) and Q = tt (open squares) subspaces. L is the 
system size. Full dots corresponds to the XY groundstate. 
The transition point is indicated in the figure. 



FIG. 2. < rik > distribution for different values of the pa- 
rameter W is plotted in (a). In (b) < A\Ak > vs k is plotted 
for differents values of W (full symbols). The open symbols 
corresponds to W — —00 (open triangles) and W — +00 
(open squares). 



FIG. 3. Zq as a function of W for Q — (a) and Q = n 
(b). The transition to jy-pairmg state occurs at W c — —1.754. 



FIG. 4. F(w) (see the text) as a function of W: It is cal- 
culated at half filling for L = 6, 8 and 10 (a). The effect of 
the band-filling for an L = 10 system with different particle 
number N e (b). 



FIG. 5. Lowest energy vs. W calculated in the Q = (open 
circles) and Q = n (open squares) subspaces. L is the system 
size and N e is the particle number. 
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